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Abstract 

Resonant photopion production with the cosmic microwave background predicts a suppression of extragalactic protons above 
the famous Greisen-Zatsepin-Kuzmin cutoff at about -Egzk ~ 5 x 10 10 GeV. Current cosmic ray data measured by the AGAS A 
and HiRes Collaborations do not unambiguously confirm the GZK cutoff and leave a window for speculations about the origin 
and chemical composition of the highest energy cosmic rays. In this work we analyze the possibility of strongly interacting 
neutrino primaries and derive model-independent quantitative requirements on the neutrino-nucleon inelastic cross section for 
a viable explanation of the cosmic ray data. Search results on weakly interacting cosmic particles from the AGASA and RICE 
experiments are taken into account simultaneously. Using a flexible parameterization of the inelastic neutrino-nucleon cross 
section we find that a combined fit of the data does not favor the Standard Model neutrino-nucleon inelastic cross section, but 
requires, at 90% confidence level, a steep increase within one energy decade around -Egzk by four orders of magnitude. We 
illustrate such an enhancement within some extensions of the Standard Model. The impact of new cosmic ray data or cosmic 
neutrino search results on this scenario, notably from the Pierre Auger Observatory soon, can be immediately evaluated within 
our approach. 



1 Introduction 

In the mid-1960's Greisen Q], Zatsepin and Kuzmin 
realized that the space-filling molasses of photons con- 
stituting the cosmic microwave background (CMB) 
limits the observation of high-energy charged particles 
originating at astrophysical sources. The attenuation 
length of protons within the CMB drops below 50 Mpc 
above the photo-pion production threshold of about 
Eqzk ~ 5 x 10 10 GeV. In the case of heavier nuclei, 
photo-disintegration with CMB photons predicts a sim- 
ilar and even stronger attenuation above this energy 
(e.g., 3]). Hence, the apparent horizon of charged ultra 
high energy (UHE) cosmic rays (CRs) is of size com- 
parable with the diameter of our local supercluster. 
Accordingly, the contribution of extragalactic charged 
particles to the CR spectrum measured on Earth is 
expected to show a cutoff at -Egzk ■ 

The UHE CR spectrum measured by the AGASA 
and HiRes [HEIGHT Collaborations is shown in Fig. 1. 
The reported flux estimates derived from UHE events 
show a qualitative difference between the data sets. 
While AGASA reports an excess of eleven events above 
10 11 GeV, the HiRes data seem to support the exis- 



tence of the GZK-cutoff. However, the significance of 
the difference is small due to the low statistics as well 
as large systematic errors in energy calibration. Apart 
from this, the origin and chemical composition of UHE 
CR remains a puzzle in astroparticle physics. The Pierre 
Auger Observatory (PAO), which will publish its first 
results in summer 2005, is expected to shed light on 
these questions. 

A recent analysis of the HiRes data 9] indicates a change 
in the UHE CR composition at around 10 8 6 GeV from 
heavy nuclei to a light component, which can be inter- 
preted as the onset of extragalactic proton dominance 
in the data. If this is correct the UHE CR data is ex- 
pected to show a cutoff at -Egzk due to resonant photo- 
pion production in the CMB. This is shown as the thin 
shaded band in Fig. 1 for our model of extragalactic pro- 
tons described in §2. On the other hand, this mechanism 
leads to the generation of high energy neutrinos from the 
decay of charged pions. This flux of cosmogenic neutri- 
nos [lOLll 1| may exceed the proton flux above and around 
-Egzk depending on the particular model parameters. In 
general, these cosmogenic neutrinos are accompanied by 
neutrinos originating from photo-hadronic processes di- 
rectly in the CR sources. Recently, two of us have been 
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Fig. 1. Top panel: The quasi- vertical flux of UHE CRs mea- 
sured by AGASA, Fly's Eye and HiRes-I/II. Bottom panel: 
The shaded band shows the flux of extragalactic protons cor- 
responding to the 90% CL of our fit. For the normalization 
we averaged the data between 10 8 ' 6 GeV and 10 11 GeV ac- 
cording to the experimental exposure (see §2). The hatched 
band shows the combined flux of cosmogenic neutrinos and 
neutrinos from optically thin sources. 

engaged in a derivation of a lower bound on this addi- 
tional source of neutrinos directly related to the observed 
flux of CRs • The sum of these two contributions is 
shown as the hatched band in Fig. 1. 

The observation that the UHE cosmic neutrino flux 
may exceed the proton flux above Eqzk has motivated 
Berezinsky and Zatsepin [lfj to speculate on the pos- 
sibility that neutrinos constitute the highest energy 
CR events, assuming them to interact so strongly as 
hadrons at the relevant energies [l^. Quantitatively, 
this requires a (rapid) rise of the neutrino-nucleon to- 
tal inelastic cross section a™ N by at least five orders 
of magnitude, above the Standard Model (SM) pre- 
diction, a™ « 7.84 pb (£„/GeV) ' 363 [H E3- The 
realization of such a behavior has been proposed abun- 
dantly in scenarios beyond the (perturbative) SM: e.g. 
arising t hrough com posi teness [l8l lip. 12(31, t hrough elec- 
troweak sphalerons [M |H llal^lalE H ill HJ, 
through string excitations in theories with a low 
string and unification scale [3(], El H3- through 
Kaluza-Klein modes from compactified extra dimen- 
sions miHmmmill, or through p-brane produc- 
tion in models with warped extra dimensions 39, 40T.l41|. 



respectively (for a recent review, see Ref . 42] ) . 

Note that neutrino-nucleon inelastic cross sections 
are in general constrained by the search results on 
normal weakly-interacting UHE neutrinos. Up to 
now, UHE cosmic neutrinos have been searched 
for in the Earth atmosphere (Fly's Eye/ AG AS A), 
in the Greenland (FORTE 143') and Antarctic ice 
sheet (AMANDA g|/RICE in the sea/lake 

(BAIKAL or in the regolith of the moon 

(GLUE ^^1). For a given flux, the low if not zero 
search results so far can be turned into model- 
independent upper bounds on the neutrino-nucleon 
inelastic cross section in the energy range where 
°f at ^ °i3v < 0.05 - 0.5 mb 48] . In strongly-interacting 
neutrino scenarios, this demands that the neutrino- 
nucleon inelastic cross section should pass very rapidly 
through this intermediate energy range. This constraint 
will be further strengthened with PAO's search results 
on deeply-penetrating showers (see e.g. |48|) , comple- 
mented by IceCub e I49l . ANITA J5fll and possibly by 
EUSO H3, SalSA gjand OWL gg. 

In this work we present a general method, which was 
previously exploited in Ref. j2^, to derive model- 
independent quantitative requirements on the cross 
section in strongly interacting neutrino scenarios. It ex- 
ploits current CR data and UHE neutrino search results, 
and can easily incorporate or be applied to other data- 
releasing experiments, notably PAO. With this method 
the central question of the viability of strongly interact- 
ing neutrino scenarios can be addressed, confining the 
requirements on astrophysics (CR sources) and particle 
physics (inelastic vN cross section) simultaneously. 

The outline of the paper is as follows. In §2 we cal- 
culate the incident fluxes of extragalactic protons and 
neutrinos assuming a power-like proton injection spec- 
trum and a cosmic evolution of the source luminosity. 
The attenuation of the initial fluxes due to interactions 
with CMB photons and adiabatic energy losses are taken 
into account with the help of propagation functions. The 
number of neutrino-induced shower and cascade events 
measured at different zenith angles will depend strongly 
on the probability for neutrino scattering on target nu- 
cleons. The effect of a varying neutrino-nucleon inelas- 
tic cross section on the measurement by different type 
of experiments will be discussed in §3. In §4 we adopt 
a flexible parameterization for the rise of the inelastic 
neutrino-nucleon cross section with energy. This allows 
us to derive quantitative criteria on the viability of the 
strongly-interacting neutrino scenarios. In §5 we present 
our conclusions and give an outlook. 

2 Extragalactic Particle Production 

Following typical models of cosmic particle accelerators, 
such as active galactic nuclei or gamma ray bursts (see 
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Fig. 2. Two possible production mechanisms of CRs by beam 
protons in a cosmic accelerator. The relative fluxes depend 
on the optical thickness of the magnetic confinement region. 

e.g. |54j for a nice review), we assume that the origin of 
extragalactic protons are /3-decayed neutrons, which are 
produced by photo-hadronic processes of beam protons 
and subsequently escape from the magnetically confined 
acceleration region. Fig. 2 shows a sketch of this mech- 
anism. The relative output of other neutral particles, in 
particular high energy neutrinos, depend on the details 
of the source such as the densities of the target photons 
and the ambient gas [H3 • In the following we will not be 
concerned about the detailed accelerator specifications 
and start directly with the injection spectrum of neu- 
trons. 

2.1 Extragalactic Protons 

The apparent isotropy of the CR data suggests the as- 
sumption of a spatially homogeneous and isotropic dis- 
tribution of the extragalactic neutron sources. Further- 
more we assume that the emissivity distribution per co- 
moving volume C n factorizes into the red-shift evolu- 
tion of the source luminosity and the injection spectrum. 
For simplicity we choose a widely used power-law ansatz 
with an exponential cut-off (see, e.g. Refs. [EH l57l l58l 
HSU): 

oc {l + z) n E-^e- E -' E — (1) 

^min ^ Z <C 2 max 

We exclude nearby (redshift z m i n ) and early (z max ) 
sources and fix these parameters in the following at 
Zmin = 0.012 (corresponding to r min sa 50Mpc) and 
-"max = 2. The maximal injection energy E max is fixed 
at 10 12 GeV in our analysis. 

The relevant flux of protons propagating towards the 
Earth resulting from the /3-decay of the source neutrons 
is attenuated due to energy red-shift, e + e~ pair pro- 
duction and pho to-pion production with CMB photons 
(see e.g. [57i l58l l59jV The last of these processes pro- 
vides the source of cosmogenic neutrinos. If we define 
P p \ n (E; E n , r) as the expected number of protons above 
an energy E given a neutron with injection energy E n 



at a distance r, then the expected proton flux incident 
on Earth can be expressed as 

oc oo 

;m 1 f at, f a 9P p \ n (E;E n ,r) 
M E ) = J dE n / dr — Cn . (2) 



The propagation functions P a \b |6ll l28j| have been cal- 
culated using the SOPHIA Monte-Carlo program 
and are available at www.desy.de/~uhecr/. The prop- 
agation distance r and the redshift z are related by 
dz = (1 + z)H(z) dr, where the Hubble expansion rate 
at a redshift z is related to the present one Hq through 
H 2 (z) = H§ \n M (l + z) 3 + n A ]. Following recent re- 
sults in cosmology 63] we will assume a flat A-dominated 
universe with relative energy densities 17m = 0.3 and 
fl\ = 0.7 as well as Hq — 72kms~ 1 Mpc -1 , the center 
value of the Hubble Space Telescope Key project 

2.2 Cosmogenic Neutrinos and Neutrinos from Opti- 
cally Thin Sources 

Isospin and charge conservation requires that each neu- 
tron produced in photohadronic interactions of beam 
protons is accompanied by a charged pion. It decays and 
produces v^, and u e 's, as sketched in Fig. 2. On aver- 
age, each neutrino carries about a quarter of the pion's 
energy. If we define €„• as the ratio E w + /E n (see Ref. 65]) 
we can express the average energy of a single neutrino 
as E v w e^En/A. In the following we will use w 0.28, 
suitable for resonant photoproduction at the energies in 
question (see Ref. [T^l. 

In optically thin sources the relative normalization of the 
neutron and neutrino emissivity is fixed by relating their 
bolometric fluxes per co-moving volume and is given by 

^C l/ (z,^E n ) = 3C n {z,E n ) . (3) 

This flux of neutrinos is directly associated with the 
injected neutrons and serve as a minimal contribution 
from the source. In optically thicker sources neutrons 
may undergo photohadronic interactions before escap- 
ing the confinement region, as sketched in Fig. 2 (right). 
This will decrease the emissivity of neutrons compared 
to that of neutrinos. Depending on the ambient gas, pp 
interactions may also dominate over photohadronic pro- 
cesses in the source and produce additional neutrinos. 

The opacity of the CMB to UHE protons propagating 
over cosmological distances guarantees a cosmogenic flux 
of neutrinos, originated in the reaction p + 7cmb — * 
N + 7T+ — > — > v tl i> ll i> e e + ... [To|, Recently, two 
of us were involved in an investigation of the actual size 
of the cosmogenic neutrino flux |5Sj . In this work, we 
consider the combined flux of cosmogenic neutrinos and 
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Fig. 3. Relative exposure 1Z (Eq. (5)) to strongly-interacting 
neutrinos as a function of the neutrino-nucleon inelastic cross 
section (t™ n , of AGASA (solid) and HiRes (dashed), respec- 
tively) . The hatched region shows the pre dicted contribution 
from SM charged current interaction |16Ul7j|. 

neutrinos from jry interactions in optically thin sources, 
which can be calculated by 

oo oo 

= ^ L / djE > / dr dE Cl ' (4) 

o o 

where the propagation functions P v \ n and P v \ v are de- 
fined analogously to P p \ n in the previous section. 

As it was discussed in Ref. this flux of high energy 
neutrinos from optically thin sources is almost in reach 
of the AMANDA-II detector for neutrino energies of the 
order 10 7 GeV and should soon be observable by its suc- 
cessor experiment IceCube. This observation is crucial 
for our assumption about the sources of CRs and the as- 
sociated flux of neutrinos from the same sources. If even 
cosmogenic neutrinos are not detected in future exper- 
iments, also the assumption on the extragalactic origin 
of CRs has to be reconsidered. 

3 Footprints of Strongly Interacting Neutrinos 

The search for UHE cosmic neutrinos is a challeng- 
ing task due to their feeble interactions in the SM. To 
overcome this problem large-scale detectors and novel 
techniques are deployed and proposed. Up to now, 
experiments with UHE neutrino sensitivity include 
AGASA, RICE, GLUE, FORTE, BAIKAL, Fly's Eye, 
AMANDA, and PAO. Further ambitious projects in- 
clude ANITA and IceCube, possibly followed by EUSO, 
SalSA, or OWL. However, if the neutrino-nucleon in- 
elastic cross section <j™ N (E) increases more rapidly with 
energy than the SM predicts, UHE cosmic neutrinos 
could already leave their footprints in the CR observa- 
tional data. This happens when the neutrino interaction 




10 9 10 10 10 11 
E [GeV] 

Fig. 4. The sensitivity in terms of the the maximal average 
flux J max per bin log 10 AE/E = ±0.05, with average cross 
section per bin a consistent to the 95% C.L. with the exper- 
imental results on QHS events at AGASA (top), contained 
events at RICE (center) and QHS events at PAO (bottom). 
The flux is shown as E 2 J max relative to the energy density 
Wcas = 8.5 x 10 6 eVm~ 2 s" 1 sr" 1 (dashed), corresponding to 
the cascade limit. The contours show the increase of the sen- 
sitivity by one order of magnitude (comp. the Appendices 
B, C and D). The dotted line is the extrapolation of SM 
charged and neutral current cross sections. 

length \ v = m p /a™ N becomes comparable with the at- 
mospheric depth x atm (0) in the quasi- vertical direction, 
e.g. < 45° (cf. Appendix A). 

In this work we aim to investigate the possibility of 
strongly interacting extragalactic neutrinos appearing as 
the highest energy CR events in dependence of a varying 
a^ N {E). We use CR data above 10 86 GeV measured by 
the AGASA and HiRes collaborations, and consider pro- 
tons and neutrino primaries only. For simplicity we as- 
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sume that the characteristics of neutrino-induced show- 
ers are indistinguishable from those induced by protons. 
In particular, we assume for both primaries i) a complete 
conversion of the incident energy into the shower, and ii) 
equal detection efficiencies at the highest energies |67| . 

Under these assumptions we define lZ(a™ N ) as the rela- 
tive experimental exposure to strongly interacting neu- 
trinos compared to the exposure £{E) to protons. The 
number of detected events is then given by 



Table 1 

Best fit (81% CL), resolution of the fit, and range of param- 
eters^ 





7 


n 


log 10 A 




lo SlO TUaV 


best fit 


2.4 


3.8 


7.0 


0.15 


10.95 


grid size 


0.01 


0.05 


0.16 


0.044 


0.1 


min 


2.00 


0.00 


0.0 


0.1 


9.90 


max 


2.99 


4.95 


7.0 


1.2 


11.50 



AUs= / dE£(E)(j p (E)+TZ(a^ N (E))J u (E)) . 

(5) 

The relative exposure 72, (o^-^)) is determined by the 
search criterion on the zenith angle 6 and the (observed) 
atmospheric depth (cf. Appendix A and Fig. A.l) 
adopted by each experiments. HiRes (9 < 60°) thus has 
a larger relative exposure than AGASA (8 < 45°), as 
Fig. 3 shows. For both experiments one sees that neu- 
trinos start to contribute significantly to quasi-vertical 
showers for a™ N > 0(1) mb. On the other hand, in the 
intermediate range er^ < a™ N < 1 mb, neutrinos can 
still penetrate deeply in the quasi-horizontal direction, 
or in the ice upper surface. For a given neutrino flux, 
the non-observation of such events so far can be turned 
into model-independent upper bounds on the neutrino- 
nucleon inelastic cross section in the intermediate range, 
or be used to constrain models which predict an anoma- 
lously enhancement of a™ N at high energies • I n the 
following we will focus on the search results on quasi- 
horizontal showers (QHS) at AGASA 68] and contained 
events at RICE [slj. Fig. 4 shows the sensitivity of 
these experiments in terms of the maximal average flux 
(J™ ax ) per bin and mean inelastic cross section {u™ N ) 
consistent with the experimental results. For details of 
the calculation see Appendices B, C and D. 

4 High Energy Inelastic Neutrino-Nucleon 
Cross Section 

For a satisfying explanation of post-GZK events by ex- 
tragalactic neutrinos, the neutrino-nucleon cross section 
has to reach nucleonic values around Eqzk as shown in 
Fig. 3. In order to be consistent with the data on QHSs 
at AGASA and events at RICE one expects only steeply 
increasing inelastic interactions as good candidates for 
a combined statistical analysis. In the following we want 
to support these rather qualitative arguments by a sta- 
tistical analysis of the existing data assuming strongly 
interacting extragalactic neutrinos. 

In frequentists statistic the level of agreement of a par- 
ticular hypothesis TL with the experimental data can be 
represented by 



the integrated probability of those samples N' which 
have a smaller probability P than the actual experimen- 
tal result A exp . In general, Tt is then accepted (or re- 
jected) at a chosen significance level Q corresponding to 
a confidence level 1 — Q [70] . As it is standard in statis- 
tics, we choose 90%, 95% and 99% as benchmarks for 
the acceptance of our model. 

In our case, the probability P is made up by Poisson dis- 
tributions of vertical events (AGASA, Fly's Eye Stereo 
and Hires-I/II; see Appendix A), QHSs at AGASA (see 
Appendix B) and contained events at RICE (see Ap- 
pendix C) with an expectation value determined by the 
hypothesis H. We also account in P for the systematic 
error in energy calibration of about 30% by a Gaussian 
distribution for a shift of the observed spectra by a mul- 
tiple of the smallest bin-size. The expectation value for 
the Poisson-distributed events is determined by the hy- 
pothesis 7i(7, n, (J™ N ), i.e. by the particular model for 
the inelastic neutrino-nucleon cross section a™ N together 
with the source luminosity given by an injection index 
7 and an evolution index n. 

The absolute value of the predicted flux is a priori un- 
known due to our lack of knowledge of the CR source lu- 
minosity. For each experiment individually we normalize 
the events induced by protons and neutrinos to the data 
between 10 8 6 GeV and 10 12 GeV. The resulting ambi- 
guity in the normalization of the proton and neutrino 
fluxes has to be removed for a prediction of horizontal 
events at AGASA and contained events at RICE. In this 
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Fig. 5. The 90% (shaded), 95% (fine-hatched) and 99% 
(coarse-hatched) CLs of the injection spectrum marginalize 
w.r.t. the neutrino-nucleon cross section. 
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Fig. 6. The 90% (shaded), 95% (fine-hatched) and 99% 
(coarse-hatched) CLs of the vN cross section. We marginal- 
ize w.r.t. the other parameters. 

case, we normalize the fluxes to an average data set in- 
terpolating between AGASA, Fly's Eye and HiRes-1/11 
according to the exposure of the individual bins. 

Now that we have set up our statistical analysis, we need 
to specify the particular neutrino-nucleon cross section 
for our hypothesis. There are various theoretical ideas 
for a rapid increase of the neutrino-nucleon cross sec- 
tion cr„7v referring to physics beyond the perturbative 
SM as was already stated in the introduction. Based on 
our previous considerations, we are interested in three 
characteristics of a strong neutrino-nucleon interaction: 



CO 




1 


10 3 


CO 




N 

1 




E 




V 








[Ge 


10 2 














Eh 





™1 

detected vs. predicted 
events at AGASA 




Fig. 7. The range of predicted events for AGASA induced 
by extragalactic protons and neutrinos corresponding to the 
90% CL of the goodness-of-fit test (shaded band) compared 
to the AGASA data shifted in energy by ±30%. 
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Fig. 8. The range of the cross section within the 99%, 
95% and 90% CL. The lines are theoretical predictions of 
an enhancement of the neutrino-nucleon cross-section by 
electroweak sphalerons [29| (short-dashed), p-branes |41| 
(long-dashed) and string excitations [3^1 (dotted). 

the threshold energy E t h marking the changeover from 
weak to strong interaction, the range zLAE/E^ for this 
change and the final amplification A of the SM predic- 
tion (Ism of charged and neutral current interactions. We 
parameterize the neutrino-nucleon cross section as 



log 



in 



>N 



Acr SM 



1 + tanh 



/ HE V /E, 



\\n(AE/E, 



th) 



(7) 



The results of our goodness-of-fit test for the inelastic 
cross section (7) are shown in Figs. 5 and 6, with the best 
fit given in Table 1. Each plot of Fig. 5 and 6 shows the 
goodness-of-fit w.r.t. two parameters as confidence lev- 
els (CLs) oil — Q. The remaining parameters of the set 
{7, n, A, AE/Eth, Eth} are marginalized by a % 2 mini- 
mization using simulated annealing as in Ref. [711. 

Comparing the boundaries of the parameters in Table 1 
with the contours in Fig. 6 it seems that the best fit is 
correlated to our choice of the maximal amplification A 
and the minimal width AE/E t h - For our statistical anal- 
ysis we kept the relative width AE/E t h always larger 
than 2 in order to account for the relative width of the 
bins of > 10 01 . The relative exposure of strongly inter- 
acting neutrinos shown in Fig. 3 indicates that the con- 
tribution of neutrinos to the vertical spectrum saturates 
under an amplification of the inelastic cross section by 
more than seven orders of magnitude. This motivates us 
to limit the amplification below A < 10 7 in our analysis. 

The goodness-of-fit test of the combined data requires, to 
the 90% CL, a steep increase by an amplification factor 
of A > 10 4 over a tiny energy interval AE/E^ < 5 
at around 10 11 GeV. In particular, the SM neutrino- 
nucleon inelastic cross section corresponding to A = 1 is 
not favored by the data. 
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Fig. 7 shows the predicted events at AGASA as a shaded 
band corresponding to the 90% CL of the fit. The ob- 
served events at AGASA are shown together with a 30% 
shift of the energy to higher and lower values. The dif- 
ference between the normalization of the AGASA and 
HiRes data can be remove d by a relative re-calibration 
of the energy scale by 30% [T^, UM ■ ^ n our approach we 
keep the normalization as stated from the AGASA and 
HiRes Collaborations and integrate the systematic error 
in energy calibration into Eq. (6) . 

For the 90%, 95% and 99% allowed range of the parame- 
ter shown in Fig. 6 we plot the range of the corresponding 
cross section in Fig. 8, which can be used as a benchmark 
test for scenarios proposing strongly interacting neutri- 
nos as a solution to the GZK puzzle. As an illustration, 
we have considered three models of a rapidly increas- 
ing neutrino-nucleon cross section based on electroweak 
sphalerons [2!| , p-branes [4l| and string excitations jlj • 

• Electroweak sphalerons : We have used the neutrino- 
nucleon cross section induced by electroweak instan- 
tons shown in Ref . [23 , based on the neutrino-parton 
cross section from Ref. 1271 . the latter exploiting nu- 
merical results from Ref. [26j . A direct fit of 7 and n 
with this cross section gives 1 — Q = 98%, which is in 
very good agreement with Fig. 8. 

• p-branes : We have calculated o~(vN — > brane) given 
in Ref. fTl] as Eqs. (9) and (10) for m = 6 extra spa- 
tial dimensions, a fundamental scale of gravity Mo = 
300 TeV and a ratio L/L* = 0.005 of small to large 
compactification radii. We let all partons interact uni- 
versally with the neutrino and use the CTEQ [74| Set 
5D parton distribution functions (PDFs) contained in 
the FORTRAN library PDFLIB [7| Version 8.04 at 
the factorization scale /i = 10 TeV. Our fit of 7 and n 
gives 1 - g = 83%. 

• String excitations : For the neutrino-quark cross sec- 
tion given in Ref. 32] with a string scale M* = 70 TeV 
and for the set of parameters No = C = 16, char- 
acterizing the width and the absolute normalization, 
respectively (cf. their Fig. 2), we derived the neutrino- 
nucleon cross section with the PDFs described in the 
previous item. Our fit gives 1 — Q = 84% for this cross 
section, again in very good agreement with Fig. 8. 

We should also mention at this point that a neutrino- 
nucleon cross section much larger than the SM predic- 
tions at some high energy scale will also have impact on 
the elastic scattering amplitude at much lower energies 
due to dispersion relations [7(j. Eq. (E.2) in Appendix E 
gives the expected relative rise of the real amplitude 

Table 2 

Relative contribution to the real part of the scattering am- 
plitude at E v — 100 GeV predicted by dispersion relations. 



best fit 


90% CL 


95% CL 


99% CL 


0.091 


< 0.13 


< 0.15 


< 0.21 











compared to the SM prediction for arbitrary models 
of the neutrino-nucleon cross section. We have checked 
that for the models displayed in Fig. 8 there is no con- 
flict with low energy data on elastic neutrino-nucleon 
scattering. Table 2 shows the maximal relative contri- 
bution of strongly interacting neutrinos to the real part 
of the elastic scattering amplitude predicted by the SM 
at E v = 100 GeV. For the 99% CL the maximal contri- 
bution is 21 %. 



5 Conclusion and Outlook 

We have shown that current data on the highest en- 
ergy cosmic rays from AGASA and HiRes may be in- 
terpreted as the combined flux of extragalactic protons 
and strongly interacting extragalactic neutrinos. For the 
flux of neutrinos associated with neutrons from opti- 
cally thin sources we derived requirements on the in- 
elastic neutrino-nucleon cross section. We found, that 
a sufficiently steep increase of the cross section within 
one energy decade around Eqzk by four orders of mag- 
nitude may serve as a model also consistent with the 
search results on quasi-horizontal showers at AGASA 
and contained events at RICE. Our result is summa- 
rized in Fig. 8 where we show the range of the enhanced 
neutrino-nucleon cross section within the 90%, 95% and 
95% CL of the fit. We have checked that the allowed re- 
gion for the cross section is compatible with theoretical 
predictions, e.g. from electroweak sphalerons, p-branes 
and string excitations. 

Our assumption of the extragalactic origin of the ultra 
high energy cosmic rays is motivated by experimental 
composition measurements. Standard mechanisms for 
the acceleration of these particles necessitates a mini- 
mal flux of extragalactic neutrinos associated with the 
observed CRs. If protons are produced in denser or op- 
tically thicker sources the flux of neutrinos is expected 
to be increased compared to the protons and the neces- 
sity of a strong inelastic cross section will be relaxed. 
This will also be the case if we extend the cut-off of the 
neutron injection spectrum E max to values larger than 
10 12 GeV, though this seems to be hard to achieve for 
astrophysical Bottom-Up sources (see e.g. [54] for a re- 
view). These ambiguities in the flux of UHE cosmic neu- 
trinos will soon be clarified by future experiments such 
as IceCube [l2|. 

The Pierre Auger Observatory will play a crucial role 
on models of strongly interacting neutrinos. Beside the 
spectrum of vertical showers with a much better statistic 
than AGASA and HiRes, the search of quasi-horizontal 
showers will soon have a stronger sensitivity to weakly 
interacting neutrinos as Fig. 4 indicates. Within our ap- 
proach it will be easy to implement any future data, 
notably from Auger, which might finally reach large 
sensitivity on strongly interacting neutrino scenarios. 
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Here, also possible correlations with distant ast rophysi - 
cal sources can give a hint on neutrino primaries |77ll78| . 



Acknowledgments 

We would like to thank Luis Anchordoqui and Steen 
Hannestad for the nice discussions, and Zoltan Fodor 
and Sandor Katz for comments and for the support of 
our numerical analysis based on their original computer 
codes. 



A Observation of Vertical Showers at AGASA 
and HiRes 

We start with the differential rate of air showers initiated 
at the point (£, 9) by particles (cosmic ray protons, cos- 
mic neutrinos etc.) incoming with energy E and of flux 
J(E) . Here, I is the distance of this point to the detec- 
tor center measured along the shower axis, and 9 is the 
angle to the zenith at the point the shower axis hits the 
Earth's surface. We will focus in the following on showers 
with axis going through the detector. For detectors us- 
ing the Fly's Eye technique this slightly underestimates 
the rate of induced event. However, for the relative ex- 
posure 1Z of strongly interacting neutrinos this should 
effect the results only for small cross sections, and then 
by a negligible amount. 

The number of induced showers due to the inelastic in- 
teraction, the cross section for which is er ln , per unit of 
length I along the shower axis, time t, area A± perpen- 
dicular to the shower axis, shower energy and solid 
angle £1 (with dQ = sin 9 d9 d(f>)) is 

_d / d 4 A ind \ 

d£ V dt dA ± d£l dE sh J 

= e^mm ^^e-^ 1 , ( A.i) 

where m p is the proton mass, and p a ir is the air den- 
sity at the altitude h. The energy deposited in a visible 
shower is related to the incident particle energy by the 
inelasticity parameter, i.e. E s h = yE. The distribution 
of y is dependent on the scattering process. 

The number of showers observed is determined by the 
trigger efficiency P(E, 9) of the detector, and the range 
of atmospheric depths within which showers induced are 
visible to the detector. The latter, after carrying out the 
integral p^i T {l,9) di = —dx, translates into a minimal 
and maximal atmospheric depth X- (9) and x + (9) , as 
shown in Fig. A.l. After integration along the line of 




10 20 30 40 50 60 70 80 
[degree] 



Fig. A.l. The range of atmospheric depths X- < x c g < x+ 
used for the calculation of Eq. (A. 2) for quasi- vertical show- 
ers at AGASA (left-hatched) and HiRes (right-hatched) 
and quasi- horizontal showers at AGASA (shaded). For the 
quasi- vertical showers we use X- pb 0kg cm -2 . 

sight of the observer, one has 

d 4 AUs 
dt dA± dfi dE Bh 

/ T i °(S)x_(9) <r in (E) x + (9) \ 

= V{E sh ,6)J{E) ie »» -e ™p . 

(A.2) 

From this we can read off the experimental exposure 
£(E) [sm 2 sr] defined as: 

A obs = f dE£(E)J(E). (A.3) 

The AGASA and HiRes spectra shown in Fig. 1 have 
been observed in the quasi-vertical direction with 9 < 
45° and 9 < 60°, respectively. In this angular range 
both detectors are sensitive to interactions in the out- 
ermost atmospheric layers and X- (9) decreases to zero. 
The maximal atmospheric depth x+(9) depend on the 
height of the detector sites and the line of sight observed 
in the atmosphere. Also, for the case of AGASA we re- 
duce the total atmospheric depth x(9) by 500gcm~ 2 in 
order to exclude those interactions too close to the de- 
tector site in order to be triggered jT^ . Note that on av- 
erage a 10 19 eV shower develops to its maximum after 
traversing 800 gem -2 in the atmosphere. 

B Quasi-horizontal Showers at AGASA 

The AGASA Collaboration has observed one quasi- 
horizontal (9 > 60°) shower (QHS) during an op- 
eration time of 1710.5 days with an expected back- 
ground of 1.72 to.mtoM (MC statistics and systematics) 
|68j . The search criteria set the following constraints 
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on the shower maximum x max : x(9) — x max (6') < 
500gcm~ 2 and x max > 2500gcm~ 2 . On average 
the shower develops its maximum after traversing 
800 gem -2 in the atmosphere. Hence the observed at- 
mospheric depth shown in Fig. A.l varies between 
X-(6) = max(17OOgcm- 2 ,x(60 - 1300 gcm~ 2 ) and 
x+(9) = x(9). The number of observed QHS events can 
then be calculated by Eq. (A. 2). As the effective detec- 
tor area for hadronic showers we took A — 56.1 km 2 . 
The horizontal detection efficiency Vhor(E) for QHSs is 
reported to be 100% above 10 10 GeV and approximately 
zero below 10 8 GeV. In between we use a power-law 
approximation cx (log 10 (-E/GeV) — 8)", which is fit- 
ted H3 to reproduce the upper bound of 3.52 events 
(95% CL) for the observation of one QHS from charged 
current interactions reported by the AGASA Collabo- 
ration. The upper plot in Fig. 4 shows the sensitivity 
of AGASA for QHSs in terms of the maximal neu- 
trino flux E 2 J™ ax (£'y) per flavor and per energy range 
log 10 Ai;/i; = ±0.05 consistent to the 95% CL with 
the observation. 

C Contained Events at RICE 

The Radio Ice Cherenkov Experiment (RICE) at the 
South Pole has searched for electro-magnetic and 
hadronic showers based on the principle of "radio co- 
herence" . During an observation time of 3500 hours 
no candidate of an neutrino-induced event has been 
observed [69j . The expected number of events can be 
approximated as: 



effective aperture has been parametrized in Ref. [48j 
through a comparison with the geometric acceptance 
published in Ref. H3- In short, to estimate the sen- 
sitivity for PAO, the following selection criteria were 
adopted: i) 75° < 9 < 90° for the zenith angle, 
ii) A max > 2500 g/cm 2 for the shower maximum, 
which corresponds to requiring X-{9) — 1700g/cm 2 
in this work. The altitude of the PAO Southern site 
(1200 m above sea level) was also taken into ac- 
count in x + {9) = x(9). For hadronic showers with 
axis falling in the array, the effective area can then 
be parametrized as Aj_(6, E) P{E), with Aj_(9, E) sa 
cosfl • 1.475 km 2 (£/eV) ai5 \ and P{E) = 1 for 
E > 10 8 - 6 GeV, while P(E) = 0.654 log 10 (B/eV)-10.9 
below this energy. The effective aperture for all showers 
(i.e. including showers with axis not going through the 
array) is roughly 1.8 to 2.5 times larger, as shown in 
Ref. B2. 

A first model-independent investigation of the sensitiv- 
ity of PAO to anomalous neutrino interactions was per- 
formed in Ref. ■ Assuming one year of operation with 
no event observed above the expected hadronic SM back- 
ground (95% CL corresponding to 3.09 events), we esti- 
mate the prospects for PAO to strengthen the existing 
constraints imposed by AGASA and RICE. The 1-year 
projected sensitivity is shown in Fig. 4 (lowest panel). 



E Dispersion Relations 



d 3 N 



dt dQ dE 



Pice 



d?UE)o? N (E)e 



Vrff(-B) 



(C.l) 



The effective detection volume V e s(E) has been deter- 
mined by MC simulations in [8l|. We approximate the 
ice target as a cylinder with V c g(E) — Inrr 2 (E) and a 
fixed height h = 1 km. From this we can approximate the 
distance d(r, 9) a quasi-horizontal neutrino has to tra- 
verse in ice before it interacts with a nucleon inside the 
detector. Hence, the depth x(9) is composed of the at- 
mospheric depth £ a tm(#) and the depth in ice rf(r, 9)p lcc . 
The sensitivity of RICE to the neutrino flux is shown as 
the center plot in Fig. 4. 

D Quasi-Horizontal Showers at PAO 

The Pierre Auger Observatory (PAO), which is actually 
comprised of two sub-observatories, is the next large- 
scale neutrino detector in operation. The Southern site 
is currently operational and growing to its final size of 
~ 3000 km 2 . 

The rate of neutrino-induced events at the ground ar- 
rays of PAO can be calculated using Eq. (A. 2). The 



Regardless of a particular model, high energy cross sec- 
tion have to fulfill criteria relying on the analyticity and 
unitarity of the S-matrix. As was emphasized in |76j the 
total cross sections at high energies are constrained by 
low energy elastic amplitudes due to dispersion relations. 
One can relate the elastic scattering amplitude [13 $tA 
to the principal value of an integral involving the total 
(anti-) neutrino- nucleon cross section cr tot = er SM +<7 new : 



UA{E V ) - 5L4(0) 



in 



dE' 



al%{2m p E') 

E'{E'-E U ) ^ E'{E' + E V 



a^(2m p E') 



(E.l) 



We assume that our high energy cross section obeys the 
Pomerantschuk theorem, i.e. c"™ — ct pn w — * for E v — * 
oo. The elastic amplitude at E v « is dominated by 
Z-boson exchange of the order of Gf/^VZ- For E v <C 
E- <C E t h — AS we can use Eq. (E.l) to estimate the 
relative contribution of new physics at low energies as: 



MA new (E v ) V2E„ 



^A SM {E V ) 0.637ttG f 



dE'- 



E> dE' \<t sm 



(E.2) 
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